A phenomenological theory of nonphotochemical laser induced nucleation 
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Our analysis of the experimental data related to nonphotochemical laser induced nucleation in 
solutions leads to the inevitable conclusion that the phase transformation is initiated by particles 
that are metallic in nature. This conclusion appears paradoxical because the final products are 
dielectric crystals. We show that the experimental results are well accounted for by the theory of 
electric field induced nucleation of metallic particles that are elongated in the direction of the field. 
However, new physical and chemical insights are required to understand the structure of the metallic 
precursor particles and the kinetics of subsequent dielectric crystallization. 

PACS numbers: 05.70.Fh, 64.60.qj, 64.70.dg, 82.60.Nh 
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There is growing experimental support for electric field 
induced nucleation of solute particles in supersaturated 
solutions. First reported by Garetz et. al. [l|, the phe- 
nomenon referred to as nonphotochemical laser induced 
nucleation (NPLIN) has been observed with both oscil- 
lating [ll-Q and static Q fields. The term 'nonphoto- 
chemical' emphasizes that there is no light absorption; 
hence, underlying changes in electronic structure capable 
of chemical reactions are ruled out. In all of the above 
cases, the final product of nucleation was found to be 
small dielectric particles. 

A clear indication that the field is the primary phase 
change driver is the alignment of nucleated particles 
along the direction of the applied field (or laser beam 
polarization). That phenomenon has led to a type of 
'polarization switching' wherein the polymorph (crystal 
structure) of the nucleated crystal can be controlled by 
applying either linear or circular polarized light 0, [^. 
The underlying mechanism remains an open question 
with many practical implications [9] . 

Our summary of NPLIN data from the literature is pre- 
sented in Fig. [l]where the peak laser intensity (or applied 
field) and exposure time are provided at which crystal- 
lization was eventually observed in solutions at various 
supersaturation levels. For our purposes, the most im- 
portant results are that: i) the field reduces the nucle- 
ation time by 13 orders of magnitude or even larger; and 
ii) it can do so at optical frequencies (~ 10^* Hz). Other 
observations include the existence of a threshold field, 
below which nucleation did not occur, and a linear type 
correlation between the cumulative fraction of samples 
nucleated and laser intensity @, @] , as well as the solute 
supersaturation @. 

The anomalous strength of the observed field effect can 
be conveniently expressed in the terms of nucleation bar- 
rier Wo that determines the nucleus induction time, r = 
To exp (Wo/fcT). Here, tq > lO^^'^ s is the characteristic 
atomic vibration or diffusion time, k is Boltzmann's con- 
stant, and T is temperature. The observed reduction of 
T by a factor of 10"^'^ requires decreasing the nucleation 
barrier by approximately AW^(£') = Wo-W{E) - SOfcT. 
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FIG. 1: Summary of NPLIN data in aqueous solutions of: 
(a) urea [l[; (b) KCl ^S]; (c) a- and 7- glycine ^2]; (d) urea 
i; (e) L-histidine i; (f) lysozyme i; (g) KCl H; (h) 7- 
glycine @]. The ordinate is the peak applied field (unless 
labeled as the threshold field) at which nucleation occurred 
within the exposure time on the abscissa. The typical laser 
wavelength was A = 1.064 ^m, except for (e) and (f) at A = 
0.532 ^m. Samples were exposed to numerous laser pulses 
(except for single pulse exposure in Ref. @) and, typically, less 
than half of the samples crystallized when irradiated even at 
the maximum intensity. 



Earlier proposed mechanisms based on the Kerr effect 1 1 
or isotropic atomic polarizability of dielectric clusters 6 1 
provided coherent qualitative features but were shown to 
produce effects on the order of only lO^'^fcT, five orders 
of magnitude below the observations. 

Here, we propose a phenomenological theory that 
shows how NPLIN evolves through nucleation of short- 
lived, progenitor metallic particles under strong electric 
field. The metallic particles have gigantic electronic sus- 
ceptibility compared to dielectric substances and plasma 
frequency much higher than the laser frequency, thereby 
allowing the polarization to adiabatically follow the oscil- 
lating laser field. Those properties are unique to metals 
and correspond, respectively, to the above cited NPLIN 



2 



results i) and ii). Consequently, we consider the con- 
cept of metallic progenitor particles an experimental fact 
rather than an hypothesis. The question of the micro- 
scopic structure of the precursor metallic particles re- 
mains to be answered and depends on the material sys- 
tem. A paradoxical feature of our theory is that it in- 
troduces metallic particles in processes which are exper- 
imentally known to result in the formation of dielectric 
crystals. 

A relevant case of gigantic polarizability is found with 
a metal needle of length H and radius R <^ H . Under 
an electric field E, it will accumulate at its ends oppo- 
site charges of absolute value q ^ EH^e corresponding 
to the dipole moment pm ~ qH = EH^e — Eefl{H / Ry , 
where e is the dielectric permittivity of the host material 
and « HR^ is the particle volume. For comparison, 
a dielectric particle of equal volume and characteristic 
dimension fl^^^ will develop, under the same field, the 
dipole moment pd ~ EAeVl Pm, where Ae is the dif- 
ference between the dielectric permittivities of the parti- 
cle (sp) and the host (e). Hence, the electrostatic energy 
gain in nucleation, pE, is higher for needle-shaped metal 
particles by a factor of {e/ Ae){H/ R)"^ > 1. 

Regarding the dynamic characteristics, we note that 
the typical metal plasma frequency in the range of ujp ~ 
10^^ — 10^^ s~^ is much greater than both the laser fre- 
quencies used in the NPLIN work and the characteristic 
dielectric relaxation frequencies corresponding to reorga- 
nization and orientation of permanent dipoles, which are 
all below ujd ^ 10^^ s^^. Therefore, the progenitor metal 
particles will behave as good metals in the laser field, 
unlike dielectric particles that would not have time to 
polarize under the optical frequency laser field. 

The above claim of strong energy gain can be made 
more quantitative in the framework of field induced nu- 
cleation (FIN) theory [HllH. FIN is a recently devel- 
oped concept of metal phase nucleation in an insulating 
host under a strong static or oscillating [ll[ field. We 
demonstrate that FIN can account for the magnitude of 
the dramatic field effect observed in the NPLIN experi- 
ments. 

In general, a particle of volume 51 in a uniform field E 
reduces the free energy according to [l4| , 
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Here, the particle polarizability depends on the effective 
permittivity, e*, and the shape of the particle through the 
depolarizing factor n. A sphere provides n = 1/3, while 
for a prolate spheroid of radius R and height H, n — 
{R/H f[\n{2H/R) - 1]. A metallic particle with Sp 
oo leads to nAe ^ e, resulting in e* = e/n, consistent 
with the intuitive estimate presented above for the dipole 
moment of a metallic needle. 

Insight to the underlying mechanism can be gained 
by estimating the effective permittivity e* that is re- 



quired to provide the observed barrier reduction AW = 
SOfcT. From Eq. setting Fe = AW yields e* = 

8tt AW / {E"^ rt) . As a rough estimate, a particle of vol- 
ume 57 ~ 1 nm'^ and the experimental E ~ 10^ V/m gives 
the requirement e* ~ 10^. In comparison, for a dielec- 
tric sphere we have e* = 3e(£p — e)/{ep + 2e), and using 
Ep = 2.2 (for KCl) and s = 1.8 (for water at A = 1.064 
^m) [31, yields e* = 0.4; five orders of magnitude too 
low. This crude estimate reveals the unlikelihood that 
a dielectric particle could provide the necessary barrier 
suppression. 

Nucleation in the presence of an electric field is de- 
scribed by the free energy. 



F = Aa-nn + FE 



(2) 



where A is the nucleus surface area, a is the coefficient 
of surface tension, and fi is the chemical potential dif- 
ference between the two phases. For the well known 
case of spherical nuclei, A = 4ttR^, ft = {A/3)ttR^, and 
the field-dependent nucleation barrier becomes Wgph = 
max{F(i?)}, 
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with Er, 
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(3) 



where Eq is the characteristic field expressed in terms 
of the classical barrier Wq — IQira^ / {'i^^) and critical 
radius i?o — 2cr//i. Note that the critical radius in the 
field is smaller than R^, 
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In FIN, nucleation proceeds through two degrees of 
freedom by forming needle-shaped particles aligned with 
the direction of the field (or beam polarization). They 
are much more efficient at reducing the electrostatic en- 
ergy because of their larger dipole moments. The exact 
shape of the elongated nucleus is not known, but mod- 
eling with either spheroidal or cylindrical particles leads 
to differences only in numerical coefficients To']. We opt 
for the mathematically more concise form of a cylindrical 
nucleus with A = 2ttRH and Q. = ttR^H, leading to the 
free energy of Eq. ^ expressed as. 
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Here we have assumed the particle to be metallic with 
£* — e/n and used the approximation n — (R/H)^. The 
contour plot in Fig. [2] illustrates how the system can 
lower its free energy more easily by forming elongated 
particles. 

The free energy of Eq. ([5]) seems to suggest that nu- 
clei with R are the most favorable. Realistically, R 
must be greater than some minimum value determined by 



3 




FIG. 2: Contours of free energy F/Wo from Eq. ([5}; posi- 
tive and negative regions are separated by tlie zero contour. 
The contour spacing is 0.5. Nucleation of elongated particles 
along path 2 over the barrier W is more efficient than nucle- 
ation over the maximum barrier of path 1. R/Ro = a is the 
minimum physically reasonable radius. 
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FIG. 3: Comparison of the field induced barrier suppression 
from Eqs. ^ and ^ for a dielectric sphere (dash-dot), 
metallic sphere (dash), and metallic cylinder (solid). Ex- 
perimentally observed barrier reduction W/Wo = 1/4 (with 
Wo = 40fcr), shown by the horizontal line, is achieved near 
3 X 10^ V/m for the metallic cylinder; in the range of ex- 
perimental data. The region Ec < E < Ec/a^ (using 
a = 0.1) is the effective range of FIN. Nucleation in the region 
E > Ecjo? is uncertain due to the requirement of ultra-small 
nuclei. The numerical values are provided in the running text. 



extraneous requirements, such as sufficient conductivity 
to support a large dipole energy or mechanical integrity. 
Based on data for other types of systems, it was estimated 
[lo| that a reasonable minimum radius is Rmin = olRq, 
where a ~ 0.1 is a phenomenological parameter. Lacking 
more concrete information, we employ the same approx- 
imation here. The free energy in the region R < Rmin 
is substantially larger than described by Eq. ([5]), since 
the energy reducing effect of the electric field cannot be 
manifested by such thin particles. 

As a simplifying approximation we consider nucleation 
along the path R/Rq = a (see Fig. [2]); alternative paths 
that start from the origin introduce only insignificant nu- 
merical factors. Then, from Eq. ([5]) the nucleation bar- 
rier and critical aspect ratio are. 
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> 1. 



(6) 



Although the nature of the metallic precursor particles 
is unknown, with reasonable values of Wq — 1 eV, Ro — 
3 nm, and a = 0.1, Eqs. ([3]) and ([6]) imply that the 
experimentally observed barrier reduction is achieved at 
a field of E ^ 10^ V/m for dielectric spheres and E ~ 
3 X 10^ V/m for metallic cylinders; well within the range 
of NPLIN data (see Fig. [T]). In a static field, e - 100 for 
aqueous solutions and the latter value is reduced to E ^ 
3 X 10^ V/m. The field dependent nucleation barriers for 
the various scenarios are shown in Fig. [3] 

Comparing Eqs. ^ and ([6]) indicates that nucleation 
of needle-shaped particles is favored when Wcyi < Wq, 
resulting in the critical field condition E > Ec = o^^'^Eq. 
The requirement on the aspect ratio He/ Rmin 3> 1 from 



Eq. ([6]) implies the upper limit E < E^j ^/a. Taken 
together, FIN is effective in the range, 



1 < EjEc < a 



(7) 



which is clearly indicated in Fig. [3l 10^ < E < 10^ V/m 
for the numerical values mentioned above. Beyond the 
upper limit {E > Ec/a^), small nuclei with R < Rmin 
are expected. The nucleation of such small particles can 
involve other physical aspects that we do not consider 
here (cf. [iBl)- Below the lower range, spherical particles 
are more probable than cylinders but the field effect is 
negligible (i.e. minimal barrier suppression). 

From the induction time t = tq exp {Wcyi/kT) and Eq. 
([51), the threshold field is given by. 
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E, 

kT In (r/ro) ' 



(8) 



It provides verifiable predictions in terms of the threshold 
field dependence on exposure time, temperature and the 
supersaturation coefficient 7 through Wq oc 1/ In^ 7 17 1 . 
Some care must be taken in interpreting the available 
data since only the upper limits of the induction times are 
known and threshold fields were only reported in Refs. |3 
and 1^. 

The above analysis was limited to the nucleation 
stage of phase transformation. However, post-nucleation 
growth (or decay) can strongly affect the number of ex- 
perimentally observed second phase particles. In gen- 
eral, the post-nucleation processes can be rather com- 
plex, including secondary nucleation of the second phase 
particles on the precursor metallic embryos, structural 
reconstruction [18|, and subsequent particle growth by 
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accretion from the solution. The reconstruction step im- 
pUes that "nucleation is, at least, a two barrier process 
in terms of the thermodynamic potential, in which the 
first barrier necessary for cluster formation is lower than 
the main barrier necessary for the transformation of the 
already formed cluster into a stable crystalline nucleus" 
[isj : it goes beyond classical nucleation theory and was 
suggested based on empirical observations. In what fol- 
lows, we attempt a qualitative description of how the 
secondary process depends on the field and solute con- 
centration. 

As illustrated in Fig. SI a newly nucleated particle 
remains unstable upon field removal unless its size has 
grown enough (above Rq in Fig. |4]) to ensure particle 
stability (continued growth) in zero field. Therefore, the 
field needs to be maintained for a certain time, Tg, to 
let a just nucleated particle evolve into the zero-field sta- 
bility region. The particle growth rate determines both 
that time and the number of stable particles found upon 
field removal. Assuming the characteristic time of field 
exposure Tp , the condition of sufficient growth takes the 
form R{tp) > Rq. 
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FIG. 4: Sketch of the particle free energy under zero field 
and strong electric field. Ro and Re show the corresponding 
nucleation barriers, while i?oo and Reo represent the radii 
above which the particle becomes energetically favorable. The 
upward arrow shows the transition that takes place upon field 
removal, which leads to particle decay. 

When a post-nucleation stage of particle formation be- 
comes the bottleneck, then the phase transformation rate 
will not be exponential in the electric field and material 
parameters, as would be typical for nucleation processes 
19j . Here, we consider a conceivable scenario inspired by 
the data in Refs. 0, @. We assume the growth stage to be 
the bottleneck that determines the number of particles 
observed upon field removal, while the characteristic nu- 
cleation time r„ <C Tp is the shortest of all the processes. 
Thus, metal nucleation takes place with certainty during 
the time Tp of field exposure. Furthermore, we consider 
the simplest hypothesis that the probability for a particle 
to grow beyond the stability radius Rq is proportional to 



the diffusion flux / of molecules from the solute to the 
particle. The latter is given by the equation (see e.g. [20| . 
p. 431) 

/ = Airr^Ddc/dr = 47ri:i(c - cooo)(i?o - Re), (9) 

where we have implied a spherical nucleus with radius 
close to the critical radius i?o, c is the solute concen- 
tration and c — cooo = Ac is the solute oversaturation. 
Because the practical fields are much lower than i?o, it 
follows from Eq. dU that Ra - Rp ~ 2RaE^/E^, which 
yields for the number of stable particles, 

N cx E'^Ac. (10) 

Note that the dependence in Eq. (fTU]) can be shown to 
hold not only for spherical particles [assumed in Eq. ©], 
but for cylindrical particles as well. Hence, when the field 
is sufficient to induce nucleation of metallic precursors 
{E > Ec), the probability of observing crystallization 
has the dependence of Eq. (fTO|) : as observed in Refs. 
andUl 

Regarding inquiry into the nature of the metallic pro- 
genitor particles, we note that FIN may provide a means 
by which otherwise chemically unstable substances can 
persist. As an example, consider the violently reac- 
tive combination of potassium in water with an enthalpy 
of reaction, AH ^ 200 kJ/mol = 2 eV/atom. Using 
At ~ 2 eV/atom - 4 x 10^ J/m^ and cr ~ 20 eV/atom 
(for structurally different phases) yields i?o ~ 20 A and 
Wo ~ fiRl - 200 eV. Then with a = 0.1, and e - 1, 
we obtain i?c ~ 3 x 10^ V/m. Although this minimum 
required field is beyond the range of the NPLIN data, 
the point is that otherwise unexpected substances may 
become stable in the presence of a sufficiently strong elec- 
tric field. In particular, one can assume metal-like parti- 
cles of some other, more complex composition than pure 
potassium with considerably lower enthalpy of reaction. 

In summary, the most pertinent aspects of the NPLIN 
data, including the magnitude of the field effect, the po- 
larization dependence, and the threshold field are well 
described within the framework of FIN. A qualitative 
picture of the post-nucleation process leads to a phase 
transformation rate that depends quadratically on the 
field and linearly on solute concentration. The paradoxi- 
cal requirement is that metallic progenitor particles pre- 
cede the formation of dielectric crystals. The structure 
of the metallic precursors and the kinetics of formation 
of the dielectric crystals are intriguing and largely open 
questions. 
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